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ABSTRACT 

We argue that both velocity-resolved reverberation mapping and spectropolarimetry require broad 
line regions (BLRs) of active galactic nuclei (AGNs) to have an inflow velocity on the order of a 
quarter of the velocity dispersion along the line of sight. We use the STOKES Monte Carlo radiative 
transfer code to show that electron and Rayleigh scattering off the BLR and torus naturally explains 
both the blueshifted profiles of high-ionization lines compared with low-ionization lines and the rest 
frame of the host galaxy, and the relative amplitudes of blueshifts between lines from different ions. 
This resolves the long-standing conflict between the inflow implied by velocity-resolved reverberation 
mapping, and the outflow implied if the blucshifting is interpreted as the result of obscuration. The 
net inflow of the BLR can only be explained if there is a significant source of viscosity, and this is 
likely to be viscosity due to the magneto-rotational instability. Energy is generated in the BLR by this 
viscous dissipation. The mass accretion rate implied by the BLR density and inflow velocity is similar 
to the accretion rate needed to power the AGN. We suggest that the BLR and outer accretion disk are 
essentially one and the same. Modelling shows that the amount of blucshifting of the high-ionization 
lines is proportional to the mass-accretion rate. This predicts that high-accretion-rate AGNs will show 
greater blueshifts of high-ionization lines, as is observed. We point out that the scattering can lead 
to systematically too high black hole mass estimates from the C IV line for some AGNs. We note the 
similarities between NLR and BLR blueshiftings, and suggest that NLR blucshiftings could also due to 
infall and scattering. Rayleigh scattering from the BLR and torus could be more important in the UV 
than electron scattering for predominantly neutral material around AGNs. The relative importance of 
Rayleigh scattering versus electron scattering can be assessed, in principle, by comparing the profiles 
of lines of the same ion, or from the same part of the emission-line region, but at different wavelengths. 
Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galax- 
ies: quasars: emission lines — line: profiles — scattering 



1. INTRODUCTION 



The structure and kinematics of the broad-line region 
(BLR) of active galactic nuclei (AGNs) has long been a 
subject of much debate, and a very wide range of struc- 
tur es and velocities have been considered (for reviews 
see iMathews fe Capriottllll985t lOsterbrock fe Mathewi 



19861: lOsterbrockl 11991 iGaske 1 et alJ fl999T 
Sulentic. Marziani. fe Dultzin-Hacvanl 20001 ). Be- 



cause of these uncertainties, it has not been clear where 
the BLR is located, what it is doing, and hence, what 
role it pla ys in the AGN phenomen o n. Reverberation 
mapp ing (Blandford fe McKed fl98l iGaskell fe Sparkel 
1986) has enabled us to probe the structures of 
the BLR and dusty torus, and in a companion paper 
(jGaskell. Klimek. fc Nazarovaff2TM hereinafter "GKN" ) 
it is argued that the BLR and torus share a similar 
flattened toroidal structure with a high covering factor 
and self shielding. 

Varying lines of observational evidence have given con- 
flicting pictures of the velocity field in AGNs . Firs t, the 
discovery of broad absorption lines (jLvndsl Il967f ) was 
unequivocal evidence that there are outflows in AGNs. 
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Radiatively-driven outflows al so naturally explained 
BLR line profiles in AGNs (jBlumenthal fc Mathews] 
|1975| ). The outflow p icture was supported by the dis- 
covery (|Gaskelll [l982f) of the blueshifting of the high- 
ionization BLR lines with respect to the low-ionization 
lines and the rest frame of the host galaxy by ~ 600 
km s . This effect required physical separation of the 
high- and low-ionization BLR clouds, a component of ra- 
dial mot ion of th e high- ionization clouds, and an opacity 
source. iGaskelll {1982) proposed that the blueshifting 
could be explained by having the high-ionization clouds 
be radially outflowing, with obscuration in the equa- 
torial plane blocking our view of the receding clouds 
(i.e., of the redshifted side of the line profile). Al- 
though the blueshifting is usually of the order of ~ 60 
km s _1 , it can exceed 4000 km s" 1 (|Corbinl I1990D . 
The blucshifting is not only found when comparing the 
profiles of high- and low-ionization lines in individual 
AGNs, but also when using spectral principal compo- 
nent analysis (SPCA) of samples of AGNs to separate 
out line profiles into independent components such as an 
"intermediate-line re gion" (ILR) and a "very broad line 
region" (VBLR) (see iBrotherton et al.lll994l ). The mag- 
nitude of the blueshiftin g is roughly in ord er of increasing 
ionization potential (|Tvtler fc Fanl [l992). It also tends 
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to be strongest in luminous, radio-quiet AG Ns (|Corbinl 
HMl iTvtler fc FarJll99l ISulentic et a D ll99^. esp ecially 
those showing broad-absorption lines (ICorbinl I1990D. 
and in AGNs with a high accretion rate (ISulentic et all 
120001: IXu et alll2n03t ILeighlv fc Moorell2004f ). The large 
blueshifts in high accretion rate AGNs have been taken as 
an in dication of strong outflowing winds in these AGNs 
( e.g.. ILeighlv fc Mooreil200l iKomossa et al.ll2008l) . 

iGaskelll (|1988l ) pointed out that the outflowing high- 
ionization BLR scenario predicted strong velocity- 
dependent time delays in the wings of the high-ionization 
lines (the blue wing should lead the red wing by twice 
the line-continuum delay), and showed that such a sig- 
nature of outflowing winds was absent at a high con- 
fidence level in velocity-resolved reverberation mapping 
of NGC 4151. This has subsequently been found to be 
the case for all other AGNs where velocity-dependent 
time delays have been studied (Koratkar & Gaskell 
19891: ICrenshaw fc Blackwelll 119901: Koratkar fc Gaskelll 
1991allbl : iKorista et ailfl995l) . The wings of BLR lines 
either vary simultaneously (as would be expected from 
Keperian motion of the clouds or isotropic motions), 
or show slight evidence for infl ow (i.e., the re d wings 
vary first). Because of this, IGaskelll (|1988l ) argued 
that motions of BLR clouds were gravitationally dom- 
inated, and hence that they could be used for determin - 
ing black hole mas s es. IPeterson fc Wandell (|1999tl2000f ). 
lOnken fc Petersonl ([2002D . and others have also showed 
that broad line widths are consistent with the r -1 / 2 fall 
off with radius expected when motions are virialized. 

There is now little doubt that Zow-ionization BLR 
clouds (i.e., those producing Mgll and the Balmcr 
lines) are predominantly orbiting in the equatorial plane, 
since their widths show the expected correlation with 
the orientation of the rotation axis ([Wills fc Brownd 
19861). disk- like line pr o files are common (e.g ., 



Eracleous fc Halpernl [l99l IGaskell fc Sneddenl [1999), 



and the orbit al motion of c l ouds h as been detected 
(IGaskelll H9M iSergeev et all 12001 iPronik fc Sereeev! 
2006). There is thus considerable confidence that black 
hole masses can reliably be estimated from /ow-ionization 
lines. However, for the /lig/i-ionization lines, the conflict 
between the kinematics implied by the blueshifted ab- 
sorption lines and the blueshifting of the emission lines 
on the one hand, and the velocity-resolved reverberation 
mapping on the other, raises serious doubts about the 
suitability of high-ionization lines such as C IV A1549 
for estimating black hole masses. Further evidence for a 
difference in the kinematics of high- and low-ionization 
lines comes from the d ifferent velocity dependencies o f 
the physical conditions ([Snedden fc Gaskelll 120041 |2008[ ) . 
The ionizing flux received by low-ionization BLR lines 
shows the dependence on velocity one would expect as a 
result of virialization and the inverse-square law, while 
the ionizing flux received by the high-ionization clouds 
appears to be independent of velocity. 

This conflict between BLR kinematic indicators can- 
not be reconciled by assuming that the low-ionization 
lines arise in a dis k while the high-ionization lines arise 
in a wind (e.g., iCollin-Souffrin et al.l [1988), because 
the vel ocity-resolved reverberation mapping s pecifically 
shows (|Gaskelllll98l Koratkar fc Gaskelll [1989] ) that the 
/iig/i-ionization BLR gas is not outflowing. Disk-wind 



models also have a problem of explaining why line profiles 
of differ ent ions are so simi lar if they have very different 
origins (jTvtler fc Fanll 19921 ). While winds certainly exist, 
and are probably energetically important, the amount of 
mass involved is small. 

In this paper we will argue that the velocity-resolved 
reverberation mapping results are correct, and that the 
entire BLR has a net inflow. In section 2 we summarize 
the evidence for a net inflow of the BLR and estimate the 
inflow velocity. In section 3 we argue that only scattering 
is consistent with producing a blueshifting when the BLR 
is inflowing, and in section 4 we use the STOKES Monte 
Carlo radiative transfer code to show that scattering off 
an infalling medium reproduces observed blueshifted line 
profiles, and also explains the dependence of blueshift- 
ing on ionization. We consider the implications of the 
infalling BLR scenario for the energy generation mech- 
anism in AGNs in section 5. In section 6 we offer an 
explanation of why high-accretion-rate AGNs (Narrow- 
line Seyfert Is = NLSls) show a stronger blueshifting, 
and point out potential systematic effects when the C IV 
emission line is used to estimate masses of high redshift 
AGNs. 

2. INFALLING BLR GAS 

2.1. The Evidence for Infall 

Although it has generally been assumed in estimat- 
ing black hole masses that the BLR clouds are in near- 
Keplerian o rbits, or quasi- random orbits (i.e., there is no 
net inflow), [Gaskell (1988|) found that infall was favored 
and that Keplerian or random orbits were excluded at the 
97% ( single-tailed) confidence level. iKoratkar fc Gaskelll 

( 1989) similarly found that non-inflowing motion was 
excluded in Fairall 9 at the ~ 95% confidence level. 
From the in tensive monitoring 1 989 IUE monitoring of 
NGC 5548 (IClavel et al.l [1991 . [Crenshaw fc Blackwelll 
(|1990f) and iDone fc Kro lik (1996J) also favored infall of 
the CIV emitting gas in NGC 5548. From the esti- 
mated errors in the lags given by ICrenshaw fc Blackwelll 

(1990) , non-infalling motion is excluded at the 93% confi- 
dence level. The 199 3 combined HST and IUE campaign 
([Korista et alj I1995D sh owed a similar deg r ee of inflow 
(see their Table 25). As IGaskell fc Sneddenl (| 19971 ) point 
out, although the statistical significance of inflowing C IV 
for one line in any one observing campaign of any indi- 
vidual AGN is not necessarily strong, the case is very 
strong when all the AGNs are considered together. 2 The 
analyse s just mentioned are mostly for CIV A 1549, but 
IGaskelll (U98l) fo und inflow of Mgll A2798 in NGC 4151 



and iWelsh et alj (|2007l) have recently found that H/3 in 
NGC 5548 is inflowing. 

Independent evidence for i nfall comes from high- 
resolution spectropolarimetry ([Smith et al.l 12005). The 
systematic change in polarization as a function of ve- 
locity across the Balmer lines requires a net inflow of a 
scattering region somewhat exterior to the Balmer lines. 

In summary, we believe that the evidence for both the 
high- and low-ionization BLR inflowing is very strong. 

2.2. The Infall Velocity 



2 Koratkar & Gaskell (1991b) also find non-statistically signifi- 
cant C IV inflow in four additional AGNs. 
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We can estimate the ratio of net inflow veloc- 
ities to random velocities along the line of sight 
from the positions of the peaks of the red-wing/blue- 
wing cross correlation functions (CCFs). The posi- 
tions of the e xpected peaks in t he wing-wing CCFs 
are m arked in iGaskelll (|1988f ) and iKoratkar fc Gaskelll 
(11989ft. The observed peak posit ions, and the one found 
by I Crenshaw fe BlackweU (l990) are all consistent with 
the net inflow velocities b eing several tim e s sma ller than 
the non-inflow velocities. iDone fc Krolikl (|1996( ) reach a 
similar conclusion from more detailed modelling of the 
velocity-dependent delays in NGC 5548. For each ob- 
ject this suggests that the net inflow velocity of CIV is 
~ 1000 km s . In a totally independen t analysis of 
the po larization structure of Balmer lines, ISmith et al.l 
(2005) suggest an inflow velocity of the scattering region 
of 900 km s -1 . We will therefore adopt an inflow velocity 
of 1000 km s" 1 . 

3. PRODUCING A BLUESHIFT FROM INFALLING GAS 

There are two main ways of producing a blueshift of 
a line from inflowing gas. O ne is by having anisotropi c 
emission from the gas clouds (Gaskcll 1982; Wilkes 1984), 
and the ot her is by having scattering off infalling ma- 
terial (e.g., lAuer fc van Blerkoml 119721 ). Although the 
emission from BLR clouds is certainly expected to be 
anisotropic, the anisotropy is muc h greater for some 
lines t han f or others. As di s cussed bv I Wilkes fc Carswelll 
(1982) and Kallm an et alJ (|1993l ). this creates a prob- 
lem in explaining all line profiles with an inflowing 
anisotropic emitting cloud model. Lyman a has particu- 
larly strongly asymmetric emission compared with other 
lines when clouds are optically thick, which is almost 
cert ainly the case for the BL R clouds of relevance here 
(see ISnedden fc Ga skcll 2007 for evidence against a sig- 
nificant optically-thick contribution to the BLR). There 
is no evide nce that Lyman a is more asymmetric than 
other lines ([Wilkes fc CarswelllfT98l ISnedden fc Gaskelll 
2008), so it would be hard for anisotropic emission from 
infalling clouds to explain the blucshifting. 

Electron scattering has long been considered to be a 
significant source o f line broadening in A GNs in genera l 
(|Kaneko fc Ohtanil[l96l lWevmaniJll97l lMathisill970h . 
and from time to time it has been i nvoked to explain 
the l i ne profiles of individual objects (jShields fc McKeel 
H981tlLaorll2006T ). 3 Scattering regions with a ne t radia l 
motion were considered bv lAuer fc van Blerkoml (fl972l ). 
If the scattering region is outflowing, scattered pho- 
tons are redshifted, while if it is inflowing, the pho- 
tons are blueshifted. Such shifts have already been 
shown in simulations of elect ron scattering in AGNs b y 
iKallman fc Krolikl (|1986h and lFerrara fc Pietrinil (|1993l ). 

4. SCATTERING IN AN INFLOWING MEDIUM 

4.1. Spherical Scattering Shells 

We have modelled the effects of an infalling scatter- 
ing medium using the STOKES Monte Carlo radiative 
transf er code, which is described in lGoosmann fc Gaskelll 
(2007). Detailed documentation, sample input, source 

3 Note howe ver, that the appa r ent ex tended wings in the AGN 
considered by [Shields & McKcc ( 19 8l]) are act ually due to red 
continuum being too low in the Baldwin (1975) spectrum (J. M. 
Shuder - private communication). 
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Fig. 1 . — Calculated shifts of the C IV line by electron or Rayleigh 
scattering. The top black solid curve shows a Lorentzian line profile 
before scattering. The other solid curves show the blueshifting 
caused by an external spherical shell of scatterers with an inflow 
velocity of 1000 km s" 1 . The areas under all curves are the same. 
In order of increasing blueshifting and decreasing peak flux, the 
curves show the effects of r = 0.5 (red), 1.0 (green), 2.0 (blue), and 
10 (purple). The dots show the blueshifting cause by a cylindrical 
distribution of scatters with r = 20 inflowing at 1000 km s _1 . 
The half opening angle is taken to be 45 degrees, and the system 
is viewed from close to the axis (i.e., through the opening in the 
cylinder.) 

code, and executables for different computer platforms 
can be freely downloaded. 4 The optical depths, r es 
to scattering by free electrons are not expected to be 
much greate r than unity inside the h igh-ionization BLR 
of an AGN. IShields fc McKeel (fl98l e stima t e r RS < 1 
for electrons between BLR clouds, and lLaorl (|2006f ) es- 
timates r es « 0.3 for typical BLR clouds. Modelling 
of BLR clouds with the photoionization code CLOUDY 
(|Ferland et aljl998j) gave similar values of r es . We there- 
fore investigated quasi-spherical external scattering re- 
gions with r = 0.5, 1, and 2. We also modelled r = 10 to 
investigate effects of significantly larger optical depths. 
The infall velocity was taken to be 1000 km s _1 in all 
models. Since in this paper we are interested in model- 
ing just the infall and line asymmetry, we simply assumed 
that the unscattered line had an intrinsic broadening due 
to bulk motions. We considered infalling scattering shells 
both inside and outside the BLR. The results are shown 
in Fig. 1. We show a comparison with the blueshifted 
CIV profiles in a typical AGN in Fig. 2, and a compar- 
ison with the VBL R line profile deduced by SPCA of a 
sample of AGNs bv iBrotherton et all (|1994r ) in Fig. 3. 

Our first result was that for a line-emitting region out- 
side an infalling scattering shell there was a negligible 
effect on the line profile (differences mostly smaller than 
the plotting lines and symbols in the figures). This thus 
verifies that when there is strong radial stratification of 
the BLR (see GKN), only the innermost high-ionization 
lines are blueshifted. This natural explanation of the 
difference in blueshift with ionization is a big advantage 
of our scattering mo del over the o utflowing-wind-plus- 
obscuration model of [Gaskcll (1982) because in the lat- 
ter model one has to contrive to have the obscuration af- 
fect the outer low-ionization lines less. In our scattering 

4 www. stokes-program. info 
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Fig. 2. — The profiles of O I A1305 (narrow symmetric pro- 
file shown in red) and CIV A1549 (thick blue line) for the quasar 
PKS 0304-392. The thin black line is the blueshifted profile pro- 
duced by a spherical distribution of scatters with r = 0.5, and 
the dashed green line is the profile produced by the same distri- 
bution with r = 1. The brown dots are the profile produced by 
the same r = 20 infalling cylindrical di stribution sho wn in Fig. 1. 
PKS 0304-392 observations taken from IWilkesI p98i) . 
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Fig. 3. — The mean intermediate line region (ILR) profiles (nar- 
row red profile), and very broad li ne region (VBLR) profiles (thick 
blue line) for the ALS sample of IBrotherton et at] |l994). The 
solid line shows the shifting caused by an external spherical shell 
of scatters with r = 2 infalling at 1000 km s _1 as shown in Fig. 1. 
The dotted green line is for a similar shell with r = 1. 



model the blueshift of a line only depends on the optical 
depth and velocity of material outside where the line is 
emitted. As discussed in section 4.4 below, we can there- 
fore easily predict blueshifting as a function of ionization. 
For example, C III] A1909 will have a blueshifting about 
half that of CIV A1549, as is observ ed to be the case 
(lCorbinlll990l : ISteidel fe Sargentjfl99l . 

It can be seen in Fig. 1 that the blueshifting and 
asymmetry increase with r. For r ~ 10 (a much greater 
electron-scattering optical depth than has been consid- 
ered for the BLR), both the blueshifting and the asym- 
metry become much greater than is observed, so we can 
rule out such high optical depths along the line of sight. 

4.2. Rayleigh Scattering in the BLR and Torus 



While large electron scattering optical depths are not 
likely in the inner BLR, our modelling of BLR clouds 
with CLOUDY showed that for clouds of sufficiently 
large neutral column densities to produce significant Fe II 
emission, Rayleigh scattering was over an order of magni- 
tude more important than electron scattering in the UV. 
This is not surprising for a region where hydrogen, which 
is the main provider of electrons, is mostly neutral, since 
a typical Rayleigh cross section is larger in the UV than 
the Thomson cross section. The potential importance 
of Rayleigh scattering in modi fying the spectra o f AGN s 
has been discu ssed in detail bv lKorista fe Ferlandl (|1998h . 
and I Led (|2005l ) has considered the effect of Rayleigh scat- 
tering on the polarization of Lyman a. 

STOKES does not explicitly handle Rayleigh scatter- 
ing at present, but the angular dependence of Rayleigh 
scattering and electron scattering is identical, and the 
variation of scattering cross section across an emission 
line is unimportant, so Rayleigh scattering can be treated 
as electron scattering. The effect of Rayleigh scattering 
off an infalling medium completely covering the source 
will be to give a very large blueshift as shown in Fig. 
1 for r ^ 10. Such Rayleigh scattering would, however, 
be accompanied by substantial blueshifted low-ionization 
absorption which is never observed. We can therefore 
strongly rule out such a spherical shell which is optically 
thick to Rayleigh scattering. 

4.3. Scattering off a Flattened BLR and Torus 

Although it is unlikely that our direct line of sight 
to the BLR has a high optical depth to electron scat- 
tering, as one goes away from the black hole the BLR 
must become optically thick and merge with the torus, 
so the optical depth to both electron and Rayleigh scat- 
tering will become subst antial in the equatorial plane 
(|Korista fe Ferlandl 1 19981). The existe nce of Compton- 
reflection humps (|Pounds et al.lll990h is evidence that 
that there are substantial electron scattering optical 
depths. Reprocessing of X-rays requires r es > 1 (see 
iGoosmann et al.ll2007f ). It is also well known that half 
of all S eyfert 2 galaxies are observed to be "Compton 
thick" faisaliti, Maiolino. fe Sal vati 1999), so most, or 
all AGNs are probably opti cally thick to electron scatter- 
ing in the equatorial plane. iGaskell. Klimek. fe Nazaroval 
(2008) have argued that in a typical AGN, the BLR, like 
the torus, has a c overi ng factor of 50% or more. As 
ISmith et all (|2004fl and ISmith et all (|2005f ) point out, a 
significant scale height of the scattering region is also 
necessary in order to explain the polarization of type- 
1 AGNs. We therefore modelled an infalling scattering 
cylindrical torus with a half-opening angle of 45 degrees, 
an optical depth of 20, and an infall velocity of 1000 km 
s . The viewing position is within the half-opening an- 
gle. We show the profiles arising from such a model by 
the dotted curves in Figs. 1 and 2. 

It can be seen that the blueshifting produced by such 
a model differs insignificantly from that produced by the 
purely spherical scattering model with r ~ 1. However, 
unlike the spherical model, the shift produced in the 
torus model depends only on the infall velocity, and is 
quite insensitive to large optical depths. This is because 
in the torus case, after one or two scatterings a photon 
has a high probability of escape within the half-opening 
angle. The main difference between the quasi-spherical 
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Fig. 4. — Mean blueshifts of emission lines (see text for details) 
versus the radii, R, of maximum emission predicted in the GKN 
model. The line is a least squares fit showing log v as a function of 
log R, 

scattering case and the torus case is that the former pro- 
duces more blueward asymmetry. For most objects, such 
as the AGN shown in Fig. 2, the difference between the 
two models is negligible. The advantages of such a flat- 
tened BLR and torus are that there is no need to fine 
tune r, and high Rayleigh scattering optical depths are 
permitted. 

4.4. Blueshifting as a Function of Emission Radius 

Our scattering model predicts that the blueshift of a 
line depends on the infall velocity immediately outside 
the radius, r, the line is produced at. The radii that 
given lines are produced at are known approximately 
from reverberation mapping, and the GKN self-shielding 
photoionization model also predicts relative radii which 
agree well with the observed radii for NGC 5548. For an 
infalling BLR we expect the net infall velocity, v oc r~ ' 5 . 
In Fig. 4 we compare the relative radii, R, lines are 
emitted at in the GKN model with the mean observed 
blueshifting velocities for AGNs. The mean blueshift ve- 
locity of C III] A1909 has been taken as half the bl ueshift- 
ing o f C IV relative to Mg II (see Fig. 3 of ICorbinl 
1990ft . and other blu eshift velocities have been taken from 



Tvtler fe Far] |l992). The predicted radii (in light days) 
are scaled to NGC 5548. While th e bluesh i fts ar e for 
the large sami j les of AGNs studied bv ICorbinl (fl990h and 
iTvtler fc Fan! (fl992h . rather than for NGC 5548 itself, 
there is a similar correlation if the mean blueshifts are 
plotted against the measured lags for NGC 5548. The 
least-squares fit line in Fig. 4 gives v oc R~ 52 which is 
consistent with the expected slope of -0.50. 

5. IMPLICATIONS 

5.1. Mass Inflow Rate 

It has long been noted that, for both pure inflow and 
pure outflow, the mass transfer r ate across the BLR is 
compa rable to the accretion rate. iPadovani fe Rafanellil 
(1988) have argued, for example, that if the BLR is 
purely inflowing, then the mass inflow rate in AGNs 
correlates well with the accretion rate. Calculating the 
mass inflow rate is straight forward. We illustrate this by 
considering the C IV emitting region of the well-studied 
AGN NGC 5548. The sizes of the emitting regions are 



known in NGC 5548 from reverberation mapping and the 
continuum shape i s also relatively well known. After red- 
denin g correction (jGaskell et al.ll2004l iGaskell fcFj cukcr 
I2008D . GKN get a bolometric luminosity of 10 45 07 ergs 
s -1 during the hig h state of NGC 5548 observed by 
iKorista et alJ (|1995l ). (This is onl y slightly greater than 
the 1 Q 44 - 85 ergs s _1 estimated by IPadovani fc Rafanellil 
1988). Using the GKN continuum, an observed CIV- 
emitting radius of 8 light-days, and an electron density 
tih = 10 10 , CLOUDY photoionization models give a 
thickness of 10 137 cm for the CIV emitting region, and 
a mass of 1 solar mass. For our adopted inflow velocity 
of 1000 km s _1 the inflow time from 8 light days is ~ 6 
years. Adopting a 50% covering factor (see GKN), this 
gives a mass inflow rate of 0.08 solar masses per year. 

If we adopt a black hole mass of 10 79 solar masses 
from the average of the many NGC 5548 bla c k hole mass 
estimates given bv IVestergaard fc Peterson! ((2006;), and 
assume a standard radiative efficiency of 10%, then the 
Eddington accretion rate is 0.8 solar masses per year, and 
the accretion rate during the high state studied by GKN 
is 0.1 solar masses per year. This is comparable to the 
accretion rate we have calculated from the C IV emitting 
gas, so we can see that the BLR mass inflow can readily 
provide the mass inflow rate needed to power NGC 5548. 
Since we know the size of the C IV emitting region well 
from reverberation mapping, and the net infall velocity 
to ~ 50% or so, the main uncertainty in this calculation is 
the density, since the mass inflow rate is proportional to 
n^ 1 , because gas emits less efficiently at lower densities. 

Although wc have calculated the mass inflow rate just 
for NGC 5548, it should be noted that there is no rea- 
son to think that NGC 554 8 is unusual in this regard. 
Padovani & Rafanclli (1988) have shown that their esti- 
mated mass flow rates (calculated assuming that the en- 
tire BLR is inflowing or outflowing) are proportional to 
the accretion rate needed to produce the bolometric lumi- 
nosity for a wide variety of AGNs, including objects such 
as I Zw 1, which we now recognize as a high- accretion- 
rate narrow- line Seyfert 1. 

5.2. Producing Inflow 

As is well known, to get a net inflow of matter there 
must be outward transport of angular momentum. This 
needs viscous coupling between gas at different radii. 
This has been a long standing problem for classical accre- 
tion disk modelling, since ordinary gas viscosity is far too 
low, and it is now recognized that the necessary viscos- 
ity comes from the magn eto- rotational instability (MRI) 
(jBalbus fc Hawlevl[l99lh . Since we have inflow in the 
BLR, we have the same problem of a need for a viscosity 
to provide the necessary angular momentum transfer in 
the BLR gas. We propose that this also comes from the 
MRI. 

5.3. The Relationship Between the BLR and the 
Classical Accretion Disk 

It is now recognized that we are viewing most type-1 
AGN s within ~ 45 degrees of the axis (see lAntonuccl 
1 19931 for a review), and hence that we view the BLR 
close to face on. Even when we see disk-like line pro- 
files, the inclinations are st ill not great (see, for example, 
lEracleous fc Harpernlll994l ). Because the observed broad 
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line widths are substantial, and the BLR has to be flat- 
tened (see GKN), there must be a substantial component 
of random BLR velocity out of the plane. Such a velocity 
is also necessary to maintain the thickness of the BLR 
needed to provide the observed covering factor. 

The overall motion of the BLR therefore has to be 
as follows: the main motion is rotational but there is 
a vertical random ("turbulent") component of velocity 
(|Osterbrockl I1978H that is almost as great as the rota- 
tional velocity. Then, as argued above, there is an ad- 
ditional, slower, net inflow. This overall structure of the 
velocity field in the BLR is indistinguishable from the 
results of accretion disk magn etohydrodynamic ( MHD ) 
calculat i ons su ch as those of IHawlev fc Kro lik ( 2001). 
IGaskelll (j2008) has argued that the size of the accre- 
tion disk is such that its outer parts (those generating 
the optical and IR emission) must extend out to within 
the BLR. We therefore propose on the basis of the sim- 
ilarities of the size, physical processes needed, kinemat- 
ics, and mass inflow rates, that the BLR and the outer 
part of the accretion disk are one and the same. If, on 
the other hand, there is a separate accretion disk in 
the equatorial plane of th e BLR (e.g., as in Fig. 6 of 
iCollin-Souffrin et al.l[l980h . there is a problem because 
the BLR clouds will be colliding with it. 

The strong magnetic fields generated by the MRI solve 
two major problems. T he first is the long-standing "con- 
finement prob l em" ( see iMathews fc Capriotti|[l985l for a 
review) . iReesI (|1987j ) has shown how magnetic fields can 
confine the BLR clouds. A second problem is the survival 
problem. Clouds with a random y elocity component 
(e.g., as proposed bv lOsterbrocklll9 78T ) have a mean time 
bet ween collisions comparable to the orbital timescale 
(see lOsterbrock fc Mathews! 11986( 1 . Cloud-cloud colli- 
sions will produce very high Mach number shocks which 
will immediately destroy the colliding clouds. Strong 
magnetic fields can prevent shocks from occurring, just 
as they prevent such collisions in MHD simulations of 
accretion disks. 

6. DISCUSSION 

6.1. High- Accretion- Rate AGNs 

As discussed in the introduction, it is well established 
that high-accretion-rate AGNs (such as NLSls) show 
strong blucshifting, and this has been interpreted as 
evidence for stro ng outflowing winds in high-accretion- 
rate AGNs (e.g., iLeighlv fc Moordl2004 : IKomossa et all 
2008). However, the blueshiftings in NLSls are merely 
the extreme of the distribution for AGNs in general, so 
there is no reason to think that they have a different 
cause from the shifts in AGNs in general. We therefore 
propose that the greater blueshifts in NLS1 imply greater 
inflow rates. This result follows quite naturally from our 
STOKES modelling shown in Fig. 1. The amplitude of 
the blueshifting obviously depends linearly on the inflow 
velocity. As was shown in Fig. 1, it also depends on the 
optical depth, especially of material near or along the 
line of sight. An inflow velocity of 1000 km s _1 , as we 
have assumed above, is adequate to explain the typical 
blueshiftings, but for the most extreme examples, such 
as Q1338+41 6, where the shift is almost 5000 km s _1 
(jCorbinlll 990( 1. a higher infall velocity is needed. In Fig. 1 
the blueshifting is roughly proportional to the product of 
the scattering optical depth and the infall velocity. For 



a given column length the optical depth depends on the 
filling factor and density. The mass inflow rate is pro- 
portional to the density and inflow rate. Thus, whether 
a greater optical depth or a greater inflow velocity is 
responsible for the increased blueshifts, the blueshift is 
proportional to the mass inflow rate. 

6.2. Estimating Black Hole Masses from C TV Widths 

There is considerable interest in estimating black holes 
masses at high redshifts. Because it is difficult to mea- 
sure H/3 at high redshift, it has been suggested that the 
FWH M of C IV A1549 can be used instead of the FWHM 
of KB (jVestergaardl 120021 : IWarner. Hamann. fc Dietrich! 
[2003D . Our conclusion that the C IV-producing region is 
infalling rather than outflowing in a wind is good news 
for such endeavors. However, the FWHM of CIV needs 
to be used with caution because, not only does scatter- 
ing cause blueshifting, but it also broadens the lines (see 
Fig. 1). The width of CIV can therefore give an over- 
estimate of the virial velocity in the CIV emitting re- 
gion. Since it is the square of the line width which enters 
into the equation for the virial mass, errors in the black 
hole mass are increased. We have noted above that the 
blueshifting increases with accretion rate, and that it has 
also been found to be correlated with radio loudness and 
the presence of broad absorption lines, so there is a dan- 
ger of systematic errors when using the FWHM of C IV 
to estimate masses. It should, however, be possible to 
empirically correct for these systematic errors by inves- 
tigating the differences in masses estimated from high- 
and low-ionization lines as a function of the blucshifting. 

6.3. Blueshifts of Narrow Lines 

It has long been known 

(|Burbidge. Burbidge. fc PrendergastHl959f ) that narrow 
lines in AGNs are blueshifted, and it has been widely 
assumed that this is a c onsequence of outflow of th e 
NLR gas and dust ( e.g., iDahari fc de Robertid Il988f) . 
IKomossa et al.l (|2008l) have shown that the blueshifting 
is particularly pronounced in high-accretion-rate NLSls 
and that the blueshiftings increase with ionization 
potential. This is strikingly similar to the situation for 
the BLR blueshiftings. We therefore suggest that NLR 
blueshiftings could also be due to infall and scattering. 

6.4. Determining the Structure of the Scattering 
Region (s) 

The ease with which scattering off inflowing material 
reproduces the observed blueshifting provides gratifying 
support for our infalling BLR picture. However, it is 
disappointing in another sense: because the blueshift- 
ing is insensitive to the geometry, and because there 
is some degeneracy between the inflow rate and opti- 
cal depth, we cannot get a strong constraint on the 
geometry and kinematics because of the limited accu- 
racy with which line profiles can be measured. For- 
tunately, other observations provide constraints on the 
geometry and kinematics. Two promising ways of 
studying the distribution and kinematics of the scat- 
ters are spectropolarimetry (e.g., ISmith et alT 2004. 2005; 
iGoosmann fc G askell 2007) a nd pola rimetric reverber- 
ation mapping ( Gaskell et al.l l2008al : Shoji et al. in 
preparation). The combination of the two methods (i.e., 
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high-resolution spectropolarimetric reverberation map- 
ping) promises to be particularly powerful. 

7. CONCLUSIONS 

We have pointed out the problems with the popu- 
lar outflow/ wind explanation for the blucshifting and 
blueward asymmetry of the high-ionization lines, and 
shown that velocity-resolved reverberation mapping im- 
plies that the broad-line region is inflowing. We have 
argued that electron or Rayleigh scattering off the in- 
falling medium is responsible for the blucshifting of high- 
ionization lines, and demonstrated using the STOKES 
Monte Carlo radiative transfer code that such scattering 
readily reproduces the blucshifts and asymmetries of the 
high-ionization lines. Our model also predicts that the 
relative blueshifts for different lines in the same AGN 
should be proportional to the inverse square root of the 
radius the lines are expected to be formed at. Available 
estimates of relative blueshiftings support this. 

An infalling BLR has many important implications. 
Viscosity is required to transport angular momentum 
outwards. As with traditional accretion disks, this must 
be due to the magneto-rotational instability. As the 
BLR inflows it releases energy. The deduced mass inflow 
rates are comparable to the mass accretion rate needed 
to power AGNs. We have therefore proposed that the 
broad-line region is the main material accreting onto the 
black hole. Taken together these conclusions suggest a 
picture where the BLR is the outer part of the accretion 
disk. 

We have argued that the magnitude of the high- 



ionization blueshifting effect is proportional to the mass- 
inflow rate. The inflowing BLR picture thus naturally 
explains why high accretion rate AGNs (NLSls) show 
the largest high- ionization- line blueshifts. 

The infalling BLR picture supports use of C IV to mea- 
sure black hole masses in high-redshift AGNs, but be- 
cause scattering broadens lines as well as blueshifting 
them, caution is necessary when using C IV line widths, 
and a systematic correction could be necessary as a func- 
tion of blueshifting. 

Similarities between NLR and BLR blueshiftings sug- 
gest that NLR blueshiftings could also due to infall and 
scattering. 

Finally, as iKorista fc Ferla nd (1998) have pointed out, 
Rayleigh scattering could be important in AGNs in the 
ultraviolet, and we have proposed a simple observational 
test for the degree to which Rayleigh scattering influences 
the blueshifting. 
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